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Abstract 

We investigate the effects of the relative dielectric coefficient on ionic flows in open 
ion channels, using mathematical analysis of reasonably general Poisson-Nernst-Planck type 
models that can include the finite sizes of ions. The value of the relative dielectric coefficient is 
of course a crucial parameter for ionic behavior in general. Using the powerful theory of 
singularly perturbed problems in applied mathematics, we show that some properties of open 
channels are quite insensitive to variation in the relative dielectric coefficient, thereby 
explaining such effects seen unexpectedly in simulations. The ratio between the total number 
of one ion species and that of another ion species, and the ratio between the flux of one ion 
species and that of another ion species do not depend significantly on the relative dielectric 
coefficient. 


1 Introduction 

For open channels, permeation and selectivity are among the most important 
properties of channel functions ([16, 17, 18, 21, 22, 31, 34, 36, 37, 38, 75, 76, 85, 91, 96]). These 
macroscopic properties are outputs of nonlinear interactions of these parameters, such as 
channel structures (channel shape and spatial distribution of permanent charges), boundary 
concentrations, electric potential differences, diffusion coefficients, dielectric properties, ionic 
sizes, and many others, through microscopic laws of electrodiffusion and chemical interactions 
(see, [3, 4, 5, 6, 10, 16, 35, 46, 47, 49, 51, 65, 73, 83, 86, 88, 90, 92, 93], etc.). 

To understand these properties, the atomic structure of the channel protein is critically 
important (see, e.g., [11, 12, 13, 16, 28, 37, 92, 104, 106, 105]). It determines qualitatively 
mechanisms of selectivity and permeation, and provides basis for high resolutions of integrating 
correlations in continuum models. This approach with structural features does not though deal 
quantitatively with the experimental measurements of current at all as can be seen by the 
absence of predicted graphs of currents or contents vs. voltage or concentration. The structures 
being discussed are determined from crystals, with somewhat different locations of atoms and 
forces (or they would not crystallize), in solutions remote from the physiological solutions in 
which the channels function, and often at temperatures around 100K. Qualitative discussion is of 
great importance but it must not be confused with attempts at quantitative models. 

The quantitative treatment of selectivity and permeation differs from the usual 
qualitative treatment in several important respects. The quantitative approach based on 
continuum models treats selectivity and permeation as expressions of the current flow through 
a channel under a variety of conditions, namely different membrane potentials, different 
concentrations and compositions of ionic solutions, and different mutations. The quantities from 
computations/analyses of models can be compared directly with experimental measurements of 
current. The quantitative models are dramatically reduced and oversimplified to be sure, but they 
are precise. Such is the nature of most physical models of condensed phases. The utility of the 
models is determined not by discussion of their oversimplifications but by their ability to fit wide 
ranges of data with single sets of parameter, and to design new systems and predict their 
properties. Physicists tend to emphasize the importance of fitting data; structural biologists tend 
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to emphasize the importance of structural detail. Both approaches are needed and such 
complement and compliment each other, in our view. 

Poisson-Nernst-Planck models were introduced to deal with the correlations produced by 
the mean electric field. Previous models assumed potentials that were constant as charges and 
structures changed. The goal of PNP was to deal with the change in potential profiles enforced 
by the Poisson part of the Maxwell equations. PNP was then refined to describe correlations of 
other types with more accuracy [4, 24, 25, 26, 27, 36, 41, 42, 43, 44, 49, 50, 51, 58, 59, 60, 61, 62, 
63, 64, 65, 72, 77, 78, 79, 80, 86, 97, 98, 103] with considerable success for calcium channels 
(particularly the ryanodine receptor) and good (but incomplete) success for the DEKA sodium 
channel. 

The essential point here is that models of the PNP type we have studied capture 
correlations well enough to explain the main experimental data from three major types of 
channels cited above (see [3] for more discussions). Ionic flow through biological channels 
exhibits extremely rich dynamics ([19, 20, 37, 38, 96], etc.). The flow is multi-scale (in both space 
and time) and the general behavior depends on many parameters such as channel structures 
(channel shape and spatial distribution of permanent charges), boundary concentrations, electric 
potential differences, diffusion coefficients, dielectric properties, ionic sizes, and many others ([3, 
4, 17, 18, 21, 22]). Some of the parameters affect each other, for example, the dielectric 
coefficients and diffusion coefficients are by no means constants and they vary as the 
environment changes. It is a great challenge to understand the full behavior of ionic flows. On 
the other hand, not all the parameters play equal role for every biologically interesting quantity. 
For example, the sign or direction of the flux of an ion species depends only on the difference 
of its electrochemical potentials at the boundary points of the channel ([27]); in particular, the 
sign of the flux does not depend on the permanent charge (either its spatial occupancy and its 
density) and specifics of the dielectric coefficients; of course, the magnitude of the flux does 
depend on the permanent charge and specifics of dielectric coefficients as well as other physical 
parameters of the biological system. Another example concerns ion channels without permanent 
charge: the shape does not play a significant role in the magnitude of individual fluxes - an 
average property of the channel shape is the main determinant in this case ([48]). 

This paper is not meant to produce complete models of any channel. Rather, motivated 
by the work in [3], we will examine a particular property of selectivity in this paper with the 
quantitative approach based on Poisson-Nernst-Planck (PNP) type models. 

In this work, we are interested in the effects of the relative dielectric coefficient on ionic 
flow properties. The value of the relative dielectric coefficient is of course a crucial parameter for 
ionic behaviors in general. From the modeling point of view, the traditional Poisson-Boltzmann 
system with constant dielectric coefficients as in the Debye-Fluckel theory has been modified 
with concentration-dependent dielectric coefficients and ion-water interaction that show 
significant improvements for the qualitative property of the mean activity coefficient (see, for 
example, [99,100]). 

We will examine some properties of ionic movements and show they are not sensitive to 
changes in the relative dielectric coefficient. For example, if we consider the ionic flow through a 
biological channel (for example, sodium channel in particular the DEKA aspartate-glutamate- 
lysine-alanine channel responsible for the rising phase of nerve action potentials) involving the 
mixture of Na + CP and K + CP, the total number #(Na)ofNa + ions and the total number # 
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(K) of K + ions in the channel could depend on the value of the relative dielectric coefficient 
significantly, but the ratio # (Na)/ # K) between these two numbers does not depend 
significantly on the value of the relative dielectric coefficient under reasonable conditions. The 
latter has been observed in extensive Monte Carlo simulations (to the surprise of the authors, 
we are told, see Fig. 8-10 in [3]). Applying the powerful singular perturbation theory to Poisson- 
Nernst-Planck (PNP) type models that include finite sizes of ions, we provide analytical 
justifications when the ratio between fluxes of different ion species and the ratio between total 
numbers of particles of different ion species of the ionic mixture do not depend significantly on 
the relative dielectric coefficient, and when the ratios might depend significantly on the relative 
dielectric coefficient. 

We emphasize that the dependences of the ratios on other physical parameters such as 
radius of the channel, protein structure, etc. are important. In the future, we will examine these 
dependences which require much more detailed analysis of PNP type models than that used in 
this paper for the present purpose. 

We stress that our PNP type models can include finite sizes of ions and allow general 
distributions of permanent charges of the channel. The key condition for our results is that the 
standard dimensionless parameter s defined in (4) is small. This dimensionless parameter s 
is a version of the Debye length divided by the length I of the channel that reflects a combined 
effect of the relative dielectric coefficient s r , the dielectric constant s 0 of vacuum, and the 
characteristic concentration C 0 of the ionic mixture. The parameter s is small mostly 
because of the length scale of the channel and the characteristic concentration of the ionic 
mixture. For general electrolytes problems, the value of s could vary several orders of 
magnitudes. Our results do not apply when the value of s is large or moderate. For ion 
channel problems, the value of s would be small when the characteristic concentration C 0 

and/or the length l of channel are large. We estimate e to be about 1CT 2 or smaller for a 
few biological situations (see the end of Section 3.1 and Remark 3.1). 

We remark that when the value of the relative dielectric coefficient is varied, we assume 
other properties such as diffusion coefficients, permanent charge distribution and shape of the 
channel remain more or less unchanged. The flexibility of these parameters would require more 
extensive analysis of PNP type models with possible modifications of the models. We will 
examine these important issues in future work. 

2 Poisson-Nernst-Planck models 


2.1 Three-dimensional model 


For an ionic mixture with n ion species, PNP reads 


V • (s r (r)s 0 V CD) = -eQj,C, + Q(r)), 

s-1 

V-Jk=0, ~J k = D k (r)C k V pi k , k = l,2, — ,n 

k B T 


( 2 . 1 ) 


where re Q. Q is a three-dimensional cylindrical-like domain representing the channel, 
Q(r) is the permanent charge density, e r (r) is the relative dielectric coefficient, e 0 is the 


3 



vacuum permittivity, e is the elementary charge, k B is the Boltzmann constant, T is the 
absolute temperature; <J> is the electric potential, and, for the k th ion species, C k is the 
concentration, z k is the valence (the number of charges per particle), ju k is the 
electrochemical potential depending on the concentrations {C ; } and finite sizes of all ion 
species, J k is the flux density, and D k (r) is the diffusion coefficient. 


2.2 Quasi-one-dimensional PNP models 

Reduction of three-dimensional PNP systems (1) to quasi-one-dimensional models was 
first proposed in [82] using the fact that ion channels have narrow cross-sections relative to their 
lengths and was rigorously justified in [68] for special cases. A quasi-one-dimensional PNP model 
is 


where 


— 1 -— (£ r (X)£ 0 A(X) — <D) = -e(Y Z S C S +Q(X)), 

A(X) dX r 0 clX j^ ss 

~~ A = 0, ~J k =-±-D k (X)A(X)cX Pt ' * = 1,2,-,n 
dX k B T dX 


( 2 . 2 ) 


Xe[0j] is the coordinate along the axis of the channel, A(X) is the area of cross- 


section of the channel over the location X . Equipped with system (2), we impose the following 
boundary conditions (see, [26] for a reasoning), for k = 1,2 


0(0) = V, C k (0) = L k > 0; 0(/) = 0, k (l) = R k > 0. (2.3) 


2.3 Rescaling of the quasi-one-dimensional PNP model 

The following rescaling (see [29]) or its variations have been widely used for convenience 
of mathematical analysis. To do so, we first introduce several quantities. 

Let C 0 be a characteristic concentration of the ion solution; for example, one can take 
C (l to be the maximal value of the boundary concentrations and the absolute value of the 
permanent charge: 

C 0 =max{L l ,L 2 ,...,L n ,R 1 ,R 2 ,...,R n , sup l£)(X)l}, 

X e[0,/] 

where sup Xe[0/1 1 Q(X) I is the least upper bound of the function l<2(X)l for X e [0,/]. Also, 
let D 0 be the maximal value of all diffusion coefficients: 

D q =max{ sup D,(X), sup D 2 (X),..., sup D n (X)}, 

Xe[0,/] Xe[0,l] Ze[0,/] 

and let s r = sup Xs[o;] ^(X). 

Using these quantities, we now make a dimensionless re-scaling of the variables in system 
(2) as follows. 
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,2 _ S r S O^B^ 

2/2 /^t 

e L C ft 


y, fc(*) 


A(X) 

/ 2 


f r W = £ 4^, D k (x)= D ^, Q(x) 

£ r D 0 

^( X) = _4_0(Z), c t (*) = ^^, 7,. 


gw 

c n : 


(2.4) 


J, 


k B T 


C A 


/C 0 D 0 


In terms of the new variables, the boundary value problem (BVP) associated to equation 
(2) and the boundary conditions (3) becomes 

' 2 d ( A „ K d(p^ 


h(x) dx 


s r (x)h(x) -f = -J^Z S C S -Q(x), 

UX y 


dJ, 


0, 


1 


^=1 

D k (x)h(x)c k d/Jk 


(2.5) 


dx k B T dx 

with boundary conditions at x = 0 and x = \ 




m= V„:= 7 ^-V,c,(0) = /,:=i; 

k B T C Q 

m= 0,c,(l )=r,:=^. 


( 2 . 6 ) 


We comment that the electrochemical potential ju k can include excess components 
ti* (beyond the ideal component ju[ d ) for the ion size effect. For example, one can include any 
excess potential of the form 

1 n 

, ex X 1 

7 rp ~ kf j 

Kb 1 j =i 

with g jk 's depending on ion sizes (see, e.g., [41, 58, 59]), or some specific approximation models 
in [32, 33], 

We also want to emphasize that our results concern a particular type of properties, that 
is, how the ratio of between fluxes of different ion species and the ratio between total numbers 
of particles of different ion species of the ionic mixture depend on variation of relative dielectric 
coefficient. We do not consider the effects of other quantities in this work. Our conclusion will 
be based on some properties, not the specific formulas which are not available in general, of the 
solutions of BVP (5) and (6). 

3 Relevant properties of the singularly perturbed BVP and its 
consequences 

A mathematical consequence of the key assumption of small s is that the BVP can be 
treated as a singularly perturbed problem. A general geometric framework for analyzing 
singularly perturbed BVP of PNP type systems (5), (6) has been developed in [26, 48, 66, 67, 69] 
for classical PNP systems and in [47, 57] for PNP systems with finite ion size. We give a brief 
discussion about the framework and refer readers to above mentioned references for details. 


5 



Due to the smallness of the parameter s , the BVP (5) and (6) has two distinct scales: 
fast scale for boundary/transition layers and slow scale for bulk behaviors. Accordingly, the 
system (5) can be split into two complementary subsystems: the fast subsystem and the slow 
subsystem. These subsystems, in the s —» 0 limits, become lower dimensional systems with the 
sum of their dimensions being that of the full system. Most importantly, these limiting lower 
dimensional systems correspond to "ideal" situations of physical problems and have a better 
chance to be more or less completely understood. This is the case for classical PNP systems ([26, 
48, 66, 67, 69]) and for PNP systems with finite ion sizes ([47, 57]). Furthermore, the geometric 
singular perturbation theory provides a way that allows one to piece together the limiting 
information from the subsystems to the full system for £>0 small. 

The essential result of the geometric singular perturbation theory for the singularly 
perturbed BVP (5) and (6) is (see, e.g., [26, 66, 67]), for small s > 0, if 

(</>(x; s),c l (x; £),c 2 (x; £),..., c n (x; s), J x (e), J 2 (s),..., J n (<?)) 

is the solution, then, for each 1 < k <n , one can express c k (x;s ) as 

c k (x\s) = c k0 (x) + ac kl (x;s), (3.1) 

where c k0 (x) , the zeroth order in s approximation, is independent of e and is continuous 
at all x except at those points where the permanent charge Q(x ) jumps (discontinuity 
points of Q{x) ). Similarly, for each 1 <k <n, 

J k (£) = J k0 +aI kl (s), (3.2) 

where J k0 is a constant independent of s. 

While quantities c k0 (x) and J k0 in (1) and (2) depend on parameters such as 
boundary conditions and permanent charge distributions, they do not depend on s . The terms 
£c kl (x;s ) and £J kx {£) depend on £, in addition to other parameters, but in higher order 
terms in £. 

We comment that, in general, for example, when permanent charges and/or ion sizes are 
present, analytical formulas for the zeroth order terms c k0 (x) and J k0 are not available. In 
special cases, such as, zero or small permanent charges and/or ion sizes, c k0 (x ) and J k0 
depend on permanent charges and ion sizes and the leading orders of c k0 (x) and J k0 in 
permanent charges and ion sizes have explicit analytic formulas ([48, 57]). 

In Figure 1 (next page), the above property for the concentrations in (1) is illustrated by 
numerical simulations using Matlab-bvp4c. (Thanks to Liwei Zhang for providing the numeric 
results.) The setup is for n = 2, z x = 1, z 2 =- 1, Q(x) = 0 for x e (0,1/3) u (2/3,1) and 
Q(x) = Q = 4 for xg (1/3,2/3) (so that the permanent charge Q(x) has a jump discontinuity 
at x = l/3 and at x= 2/3), and £ = 0.05 for the profiles in the left panel and £ = 0.01 for 
the profiles in the right panel. The boundary conditions in (6) are taken to be V 0 = -0.5 , 
/, = l 2 = L - 2 , r x = r 2 = R = 3 

The figure shows numerical solutions of BVP associated to the PNP system (5) with ideal 
electrochemical potential and the boundary condition (6), where (/>{x\£) is the starred curve. 
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Cj(x;£) is the solid curve and c 2 (x;£ )) is the dashed curve. The profiles of </>(x;£), Cj(x;£) 
and c 2 (x;s ) have layers around the discontinuity points jc = 1/3 and x = 2/3 of Q(x). As 
s changes from 0.05 in the left panel to a smaller value 0.01 in the right panel, the layers 
become sharp. In the £ —>■ 0 limits, the orbits of C[(x;£) and c 2 (x;s ) in (1) would limit to 




Figure 1: A numerical solution of BVP (5) and (6) with n = 2, z x - 1, 
z 2 = ~ 1 £ = 0.05 (left panel) and £ = 0.01 (right panel), e r (x)-l , 
Q(x) = 0 for x e (0,1/3) u (2/3,1) and Q(x) = Q = 4 for xe (1/3,2/3) , 
V 0 = -0.5, /, = l 2 = L = 2 and r l =r 2 = R = 3 : (f)(x\£) is the starred curve, 
c^x',6) is the solid curve, and c 2 (x;s ) is the dashed curve. Two values of e 
are used to illustrate the trend of sharpness of the layers as s changes from 
0.05 for the left panel to a smaller value 0.01 for the right panel. 


those of c 10 (.x:) and c 2Q (x ) with vertical line segments as limiting layers; in particular, c 10 (;t) 
and c 20 (x) would satisfy the electroneutrality condition zyc M (x) +z 2 c 2Q (x) + Q(x) - 0 at every 
x except a: = 1/3 and x = 2/3 where Q(x) has a jump discontinuity, and c 10 (x) and 
c 20 (x) have jump discontinuities at, and only at, x = 1/3 and x = 2/3. 

Note that, in (1) and (2), the quantities c k (x\s) and / k (s) are scaled quantities in (4). 
In terms of the original unsealed quantities, we have 


c k (X, £ ) — CqCj, (x, £) — C Q C kQ (x) + £C 0 C ki (x, £), 
/,(£) = lC 0 D 0 J k (£) = lC 0 D 0 J k0 + £lC 0 D 0 J kl (£). 

Concerning C k ., it is noted from (4) that, when C () is large. 


£ = 


I £ r £ 0 


e 2 l 2 C r , 


is small, but 


(3.3) 
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sC n 


s r £ 0 k B TC 0 


e l l 


ill 


is not small in general (although, relative to C 0 , it is small). Therefore, ignoring the second term 
sC 0 c kl (x;s ) in (3) for the value of C k (x;s ) is tricky. Similar remark applies to the flux J k (s). 


3.1 Ratio of total numbers of different ion species 

We now apply the results in (3) to the ratio of total numbers of two different ion species. For 
definiteness, we take Na + and K + ions. 

Let c x (x\s) and c 2 {x\s) be the concentrations of Na + and K + in terms of the scaled 

variables, respectively. Then, the total numbers of Na + and K + ions in the channel are, 
respectively, 

N l (e)=tc i (X\e)dX and N 2 (e) = fc 2 (X;e)dX. 

JO JO 

It follows from X = lx and C k (X;e) = C Q c k (x;e) in (4) that 

N 1 (s)= f C 1 (X;e)dX = f lC 0 c 1 (x;e)dx, 

Jo Jo (3.4) 

N 2 (e) = | C 2 (X;s)dX = ^lC 0 c 2 (x',e)dx. 


From extensive Monte Carlo simulations for DEKA Na channel in [3], it shows (Fig. 8-10 in 
[3]) that the ratio N l /N 2 decreases as the channel diameter increases BUT the ratio does not 
depend on s r significantly, say, as s r changes from 1 to 80 . We now provide an 
explanation to the latter from formulas (1) and (3). 

It follows from (1) or (3) that 


A/ - , is) = IC 0 f c ]0 (x)dx + slC 0 f c n (x; s)dx, 

JO Jo 

N 2 (s) = 1C 0 [ c 20 (x)dx + dC 0 [ c 2l (x;e)dx, 

JO JO 

and hence, the ratio between the total number of Na + and that of K + is 

N,(e) 

N 2 (s) 


(3.5) 


ZC 0 | 

j ’ c l0 (x)dx + slC 0 

•1 

c n (x\s)dx 

‘Co J 

•1 

c 20 (x)dx + slC 0 

0 • 

•1 

c 2l (x;£)dx 


| c 10 (x)dx + s J c, | (x; s)dx 
Jc 20 (x)dx + c 2l (x; s)dx 


(3.6) 
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+ 0 ( 8 2 ). 


Note that the scaled variables c k0 and c kl are of order 0(1) quantities. If the 
denominators in (6) are not small, that is. 


[ c 20 (x)dx > 0 is not very small relativeto s, 

JO 

then the coefficient of the s order term is of order 0(1). In particular. 


N^s) 

N 2 (s) 



(3.7) 


(3.8) 


with the right-hand side independent of s and with an error bound in the order 0(s ). Note 
that, in the derivation (6), the large quantity C 0 is dropped when ratio N l (s)/N 2 (s) is 
considered. 

We now use the following specific values of the quantities involved to estimate the value 
of s in (4) for an assessment of the error bound. 

• The length between the two baths containing the channel is l(nm) = /xlO 9 (m); 

• e « 1.60 xlO -19 (C), s 0 « 8.85 xlCT 12 (F-m 1 ); 

• k B *1.38xKT 23 (JK” 1 ), T* 273 (K). 


Then, from (4), 


£ 2 = 


k B Ts r e 0 


l 2 xlO 


-18 


,3 2 


xC 0 xN A xlO e 


l-c„ 


If we take 1 = 2.5 (the length of the channel would be 2.5 (nm)) and C 0 = 10 (the 
maximum of concentrations including the concentration of the permanent charge would be 10 
(M)), then 

x5. 84xl0~ 3 (3.9) 

As relative dielectric coefficient s r changes from 1 to 80, the value of s changes by a 
factor of ^^80 . If the length / of the channel is 2.5 (nm) and the maximum of the 

concentration C 0 of the ionic solution is 10 (M), then s is of order 10~ 3 (see (9)) so a 
factor of 9 in the variation of s r makes about 5% difference in the approximation (8) and 
is not significant. 

Remark 3.1 Our quasi-one-dimensional model of a channel setup includes access regions and 
baths, following [68, 82]. The variable diameter (equivalently, the variable cross sectional area 
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A(X) in (2)) of the 'channel' in this setup includes a region with the diameter of the actual pore 
of the channel protein, including the taper often seen in structures as they join the extracellular 
space. In our setup, the diameter of the 'channel' is dramatically increased outside the channel 
protein to mimic (crudely but adequately for most purposes) the huge increase in the cross 
sectional area available for current flow, as ions leave the channel protein and enter the 
extracellular space with its macroscopic cross sectional area. What is important is that resistance 
to electric current flow and effective 'resistance' to diffusional flow is millions of times less in the 
bath region than inside the channel protein. Thus, in this treatment the quantity I is the length 
between the two electrodes where boundary conditions (3) are imposed. The channel would 
correspond to a subinterval of [0,1], or in the scaled variable x, to a subinterval, say [a,b], of 
the interval [0,1J. For this setting, cfx;£) and c 2 (x\s ) would be different from those in the 
previous setting. The total numbers of Na + and K + ions in the channel would be, respectively, 

rb rb 

N x (s) = \ lC Q cfx\s)dx and N 2 (s) = \ lC Q c 2 (x;s)dx. 


Following the lines of derivation in (6), the ratio between the total number of Na + and that of 
K + would be 


Hence, 


N 2 (s) 


rb rb rb rb rb 

c lQ (x)dx c n dx\ c 2Q dx- c w dx\ c 2l dx 

Ja | Ja Ja Ja _ Ja _hO(£’^) 


rb 

c 20 (x)dx 

Ja 




<f 


c 20 (x)dx) 


rb 

TV, (<£■) J fl c io ( x )dx 
N 2 (s) f c 20 (x)dx 

Ja 


(3.10) 


with the right-hand side independent of s and with an error bound in the order O(s). 

Since the quantity l is taken to be the length between the two electrodes where the 
boundary conditions (3) are imposed, it would be much larger than the length of the channel. 
For example, if we take l = 10 (nm), which is four times the length 2.5 (nm) used for I in 
the estimation in (9), then the quantity s in (9) would be reduced by a factor of 1/4 . As a 
result, the factor of 9 caused by the change of s r from 1 to 80 makes about 1.25% 
difference in the approximation (10). 


When conclusion (8) might fail. Recall that our conclusion, from (6) to (8), that the ratio does not 
depend significantly on the dielectric coefficient generally requires also condition (7), in addition 
to s being small. It is clearly possible that condition (7) does not hold; for example, in [3], when 
selectivity of Na + vs. Ca 2+ was examined for DEKA Na channel, a small amount (1 mM) of CaCb 
was added to NaCI (100 mM) bulk solution so that Ca-concentration inside the channel is of order 
0(10 3 ) relative to Na-concentration (see Fig. 3 in [3]). In this case, the conclusion (8) might fail. 
This is not surprise and an explanation for this possible failure goes as follows. Suppose, on the 
contrary, that the conclusion (8) still holds for this case. Since the numerator of the right-hand 
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side of (8) is of order 0(1) due to the dimensionless scaling, the denominator would be order 
0(1O“ 3 ) which is comparable to a; that is, the condition (7) does not hold. Going back to 
formula (6), the second term, due to its much smaller denominator, might not be small relative 
to the first term anymore. In this situation, the approximation (8) for the ratio fails in general; 
that is, the ratio would depend on the relative dielectric coefficient in more significant way. In 
fact, this is indeed the case as shown in Fig. 10 in [3]. This is consistent with the requirement of 
(7) in our result. Looking closer at Fig. 10 in [3], from Fig. 10 A, one sees that for a - 80, the ratio 
is not large which indicates that (7) might hold (although not always since the numerator of the 
right-hand side of (8) might be small so that the denominator could still be small). This would 
imply our conclusion holds, that is, for a near 80 , the ratio does not depend on a 
significantly. This statement agrees well with the result shown in Fig. 10 B over the region 
[40,80] for a which shows the slope of the ratio against a over this region is less than 
(4-1.5)/(80 -40) ~ 6% . 

3.2 Ratio of fluxes of different ion species 

For the same reason, the ratio of fluxes of two different ion species does not depend on 
the value of dielectric coefficient in a significant way. Indeed, if J x (a ) and J 2 (a ) arethefluxes 
of two ion species, then, from (3), 


_ lC 0 D 0 J 10 + alC 0 D 0 J u (a) _ / 1 q + £/ 11 (£') 
J 2 (a) lC 0 D 0 J 20 + alC 0 D 0 J 2X (a) J 20 + aJ 2X (a) 


10 , „ A 1^20 ^ 10^21 , /-*/ „2 > 


+ a- 


20 


J 


+ 0(a ). 


20 


Similar to condition (7), 


if I J 20 I is not small relative to a. 


(3.11) 


(3.12) 


then 

J \ ( g ) ^ "^10 
J 2 (a) J 2 0 

with the right-hand side independent of a and with an error bound in the order of 0(a). 
Hence, when a is small, the ratio of the two fluxes does not depend on the value of dielectric 
coefficient significantly. 


4 Conclusion 

Some continuum theories of ion conduction allow reasonably general conclusions 
concerning physiological parameters and their dependence on underlying physical parameters. 
The PNP theory of open ionic channels can be analyzed by quite general mathematics to 
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understand the dependence of selectivity on the relative dielectric coefficients. In this short 
paper, based on PNP type models that allow finite sizes of ions, we provide theoretical 
justifications when the ratio of the total numbers of two ion species does not depend on the 
relative dielectric coefficient in a significant way and when it might. A similar result is established 
for the ratio of fluxes of two ion species. In particular, the result explains that the selectivity does 
not depend on the relative dielectric coefficients in a significant way for ionic mixtures of K + CI" 
and Na + CI" but may depend on the relative dielectric coefficients in a significant way for ionic 
mixtures of Na + CI" and Ca ++ Cl 2 which has been observed from extensive Monte Carlo simulations 
for DEKA Na channel in [3], We believe that PNP theory of open ionic channels can provide an 
alternative way to understand ion channel properties. 
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